A Ghoshal-like Test of Equilibration in Near-Fermi-Energy Heavy Ion Collisions 
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Calorimetric and coalescence techniques have been employed to probe equilibration for hot nuclei 
produced in heavy ion collisions of 35 to 55 MeV/u projectiles with medium mass targets. Entrance 
channel mass asymmetries and energies were selected in order that very hot composite nuclei of 
similar mass and excitation would remain after early stage pre-equilibrium particle emission. Inter- 
comparison of the properties and de-excitation patterns for these different systems provides evidence 
for the production of hot nuclei with decay patterns relatively independent of the specific entrance 
channel. 
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INTRODUCTION 



Determining the extent to which the properties and 
subsequent de-excitation of an excited nucleus depend 
upon the entrance channel characteristics of the collision 
used to prepare that nucleus is a time honored way of ex- 
ploring equilibration in nuclear systems. In the early, low 
energy, experiments of Ghoshal 0, indications that the 
decay of identical compound nuclei, produced at the same 
excitation energy in reactions induced by protons and al- 
pha particles was independent of the entrance channel 
provided very strong evidence for the validity of the in- 
dependence hypothesis 0, 0], in statistical decay theo- 
ries. In higher energy collisions at near Fermi energies, 
in which faster emission processes remove energy and 
mass and the subsequent composite nucleus undergoes 
some dynamic evolution prior to a possible equilibration, 
the establishment of equilibration and entrance channel 
independence is more difficult. Nevertheless, it is pre- 
cisely in this domain that nuclei are being investigated 
to determine whether or not their disassembly modes 
are dynamically or thermodynamically driven and, in the 
latter case, whether they provide evidence of liquid-gas 
phase changes and critical behavior common to other liq- 



uids 0, H, IE 0- Tests of the degree of equilibration 
achieved in these Fermi energy collisions are clearly of 
interest. 

In this paper we report on a Ghoshal-like test of the 
equilibration of medium mass hot composite nuclei pro- 
duced in central collisions of 35 MeV/u "^^Zn-h ^^Mo, 40 
MeV/u 4°Ar -h "^Sn and 55 MeV/u^^Al + Sn and in 
mid-central collisions of 47 MeV/u ^"^Zn -t- ^^Mo. Cor- 
relations between ejectile energy and emission time pre- 
dicted by AMD-V transport model calculations 0, 0| 
were employed to follow the evolution of the tempera- 
tures from the time of first particle emission to the time 
of equilibration. The early temperature evolutions for 
the four different systems show a close correspondence 
although the time range for the observed evolution varies 
somewhat with reaction system. A very similar double 
isotope ratio temperature, Thhc ~ 6 MeV, is reached at 
the end of the pre-equilibrium stage of each reaction. The 
average masses and excitation energies remaining at the 
end of the pre-equilibrium stage were determined using 
calorimetric techniques. The derived masses range from 
85 to 114 amu and the excitation energies range from 4.88 
to 5.31 MeV/u. Comparisons of the subsequent evapo- 
ration spectra and particle multiplicities for the four sys- 
tems provide strong evidence that the de-excitation cas- 
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cades for these systems are very much ahke. The com- 
bined results indicate that very similar equilibrated hot 
nuclei are produced in the four different collisions. 

EXPERIMENT 

The reactions 35 MeV/u 64Zn+ ^^Mo, 40 MeV/u "OAr 
+ "^Sn, 55 MeV/u ^^Al + ^^'^Su and 47 MeV/u ^^Zn 
+ ^^Mo were studied at the K-500 super-conducting cy- 
clotron facility at Texas A & M University. The first 
three reaction systems were selected because various 
model calculations predict that central collisions lead, 
after pre-equilibrium emission, to hot nuclei with very 
similar composite masses and excitation energies. The 
^^Al -|- -'^^''Sn reaction has a significantly different N/Z 
ratio than the first two reactions. For the fourth reac- 
tion, 47 MeV/u ^'^Zn + ^^Mo, selection of mid-central 
collisions was predicted to lead to a composite system 
of comparable energy and mass to those of the central 
collisions for the first three systems. 

For these studies we used the NIMROD detector ar- 
ray which consists of a An charged particle array inside 
a 47r neutron calorimeter. The charged particle detec- 
tor array of NIMROD includes 166 individual Csl de- 
tectors arranged in 12 rings in polar angles from ^ 4° 
to ~ 160°. In these experiments each forward ring in- 
cluded two "super-telescopes" composed of two Si detec- 
tors and 7 Si-Csl telescopes to identify intermediate mass 
fragments (IMF). The NIMROD neutron ball, which sur- 
rounds the charged particle array, was used to determine 
the neutron multiplicities by thermalization and capture 
of emitted neutrons. As the neutron ball provides little 
neutron energy or angle information, in this experiment 
we also employed five discrete liquid scintillator neutron 
detectors of the Belgian- French DEMON array in or- 
der to obtain neutron energy and angular distributions by 
time of flight techniques. These neutron detectors were 
placed at laboratory angles of 30° to 110° relative to the 
beam direction and at distances of 2 to 3 meters from 
the target. For these measurements the top quadrant of 
the mid-section of the neutron ball was raised to allow 
5 cm gaps between neighboring sections. The DEMON 
detectors viewed the target through one of these gaps. 
Given the large amount of material in the NIMROD de- 
tector, this configuration is not ideal for measurements 
of discrete neutron spectra. However the neutron ball 
is a very effective absorber for scattered neutrons. By 
employing a standard reaction, 26 MeV/u ^°Ar -|- ^^Mo, 
for which neutron spectra and multiplicities have pre- 
viously been well established by Yoshida et al. ^3|> '^^ 
empirically determined the required DEMON detector 
efficiency curves appropriate to our experimental config- 
uration. Much greater detail on the detection systems 
and calibrations may be found in references 

During the experiment, data were taken employing two 



different trigger modes. One was a minimum bias trig- 
ger in which at least one of the Csl detectors detected a 
particle. The other was a high multiplicity trigger which 
required detected particles in 3 to 5 Csl detectors (de- 
pending upon the reaction studied). 

DATA ANALYSIS 

Many of the techniques applied in this analysis have 
been discussed pre viously in greater detail in refer- 
ences H ES H m E m. Only a brief summary of 
these techniques is included in the present work. 

An inspection of the two dimensional arrays depict- 
ing the detected correlation between charged particle 
multiplicity and neutron multiplicity in NIMROD (not 
shown) , reveals a distinct correlation in which increasing 
charged particle multiplicity is associated with increas- 
ing neutron multiplicity. Although there are significant 
fluctuations reflecting both the competition between the 
different decay modes and the detection efficiencies, these 
correlations provide a means for selecting the most vio- 
lent collisions. For the analysis reported in this paper 
we have selected events corresponding to the largest ob- 
served neutron and charged particle multiplicities. For 
35 MeV/u 64Zn-K ^^Mo, 40 MeV/u ^°Ar + "^Sn and 55 
MeV/u^^'Al + i24Sn, the 10% highest multipHcity events 
taken with the minimum bias trigger were selected in this 
manner. This selection emphasizes the lower impact pa- 
rameter collisions. For the 47 MeV/u ^^Zn + ^^Mo we 
selected a range of mid-central collisions having charged 
particle and neutron multiplicities similar to those for the 
events selected for the 35 MeV/u ^'^Zn+ ^^Mo reaction 
in order to have a comparable excitation energy. 

For the selected events we carried out analyses using 
three source fits to the observed energy and angular dis- 
tributions of the light charged particles. The assumed 
sources are the PLF ( projectile-like fragment ) source, 
the target-like fragment source (TLF) a.nd an intermedi- 
ate velocity (IV) source IHIllIllli^. As in the ear- 
lier work, the IV source typically has a source velocity 
very close to half of that of the projectile, i.e., that ex- 
pected for nucleon-nucleon scattering. From these fits we 
obtained parameters describing the ejectile spectra and 
multiplicities which can be associated to the three differ- 
ent sources. Given the continuous dynamic evolution of 
the system, such source fits should be considered as pro- 
viding only a schematic picture of the emission process. 
We have employed them to estimate the multiplicities 
and energy emission at each stage of the reaction. 

For the four reactions studied, the parameters for emis- 
sion of p, d, t, ■^He and ^He from the different sources, 
derived from the fits, follow the trends of earlier reported 
values at such projectile energies [lO, [iTl lla llH 113 ■ ^^''^ 
each individual reaction system the IV source slope tem- 
peratures for p.d,t '^He and ^He are quite similar. They 
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have values in the range of 11 to 16 MeV, characteris- 
tic of those for pr e-eq uihbrium emission in this projectile 
energy range [ifl llS llM llOl ■ The slope parameters 
for the TLF sources are much lower, in the range of 4-6 
MeV. For both sources the observed spectra result from 
a summation of the spectra of particles emitted over a 
range of time. Thus the observed slope temperature val- 
ues are affected by the relative emission probabilities over 
that time period. 



EXCITATION ENERGY DETERMINATIONS 

The masses and excitation energies of the hot nuclei 
which remain after the early (PLF and IV source) emis- 
sion have been determined from the source fit parame- 
ters. The masses were obtained by subtracting, from the 
total entrance channel mass, the mass removed by pro- 
jectile source and intermediate source particles. These 
were determined from integrations of the source fits over 
An. 

For each reaction system the excitation energy remain- 
ing in the TLF source was determined using calorimet- 
ric techniques. The data allow two different methods to 
make such a determination. The first is by subtraction, 
from the initial available energy, of the total energy (ki- 
netic energy and Q value) removed by the PLF and IV 
sources and the kinetic energy of the primary remnant. 
For this purpose, the primary remnant kinetic energy 
was determined from the TLF source velocity and the 
mass not removed by the IV and projectile like sources. 
The second method consists of a reconstruction of the 
TLF excitation energy by using the information on the 
multiplicities and kinetic energies of light charged parti- 
cles, neutrons and intermediate mass fragments associ- 
ated with the TLF source. 

E* =Y, M,pii)Ecpii) + MnEn + Q + E^ (1) 

i 

where the sum extends over charged particle type i, 
Mcp{i) and Mn are the average multiplicities of charged 
particles and neutrons, Ecpii) and i5„ are their aver- 
age kinetic energies, E^ is the total energy of the emit- 
ted gamma rays (estimated as 10 MeV) and Q is the 
Q value for the observed de-excitation starting from the 
primary TLF and assuming IMF masses consistent with 
the EFAX parameterization jij]. Some differences were 
observed in the results from the two different methods 
of excitation energy determination. These reflect inefh- 
ciencies in IMF detection resulting from incomplete cov- 
erage with Si detectors, energy thresholds and emission 
at small laboratory angles not subtended by the detec- 
tors. As the two methods should agree we employed an 
iterative analysis to achieve convergence for the two dif- 
ferent techniques, using the missing IMF multiplicity as 



a variable. For this purpose, the average masses and ki- 
netic energies of these missing IMF were estimated from 
results of earlier experiments, on some of the same sys- 
tems 

A summary of the mass and excitation energy deter- 
minations for the TLF sources in these four systems is 
presented in Table 1. The derived excitation energies 
range from 4.88 to 5.31 MeV/u with estimated overall 
uncertainties of ±10%. In the following discussion we 
take the derived excitation energies to be those of the 
equilibrated expanded systems [23, , at the time cor- 
responding to the latest particle emission associated with 
the intermediate source and thus appropriate to the tem- 
perature determined at that time. 

TEMPERATURES AT THE END OF 
PRE-EQUILIBRIUM EMISSION 

We recently reported on the use of double isotope yield 
ratio measurements to determine the temperature evolu- 
tion in intermediate energy heavy ion collisions |inj |. At 
intermediate energies the observed spectral slope param- 
eters derived from the source fits represent only some 
weighted average values as the observed spectra are con- 
volutions of the spectra at different emission times. We 
have therefore employed isotope yields to determine the 
temperatures at the end of the pre-equilibrium emission 
stage of the reaction. Using the same techniques as in ref- 
erence we have determined the double isotope yield 
ratio temperatures as a function of ejectile velocity for 
the four different systems under consideration. The ve- 
locity employed is the "surface velocity", Vsurf, of the 
emitted particles, defined as the velocity of an emitted 
species at the nuclear surface, prior to acceleration in 
the Coulomb field 0|. The energy prior to Coulomb 
acceleration is obtained in our analysis by subtraction 
of the Coulomb barrier energy derived from the source 
fits. Since the early emitted light particle energies are 
strongly correlated with emission times, and the evapo- 
rative or secondary emission contributions to the spectra 
are primarily at the lower kinetic energies, the yields of 
higher energy particles are relatively uncontaminated by 
later emission processes. 

The double isotope ratio temperatures employed are 
T^HHe, derived from the yields of d, t, ^He and ''He clus- 
ters. To calibrate the time-scale associated with our data 
we have used the results of AMD -V transport model cal- 
culations to determine the relationship between calcu- 
lated emission time and average ejected nucleon velocities 
for the different systems. The small charged ejectiles - d, 
t, '^He and ^He, emitted at early times are treated as re- 
sulting from coalescence of the nucleons 0, 0, Q, 0| . 
To focus on the earlier evolution of the system we selected 
such clusters emitted at mid-rapidity, i.e., at angles of 70 
to 80 degrees in the IV source frame. In this way we 
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attempted to isolate the emission associated with the IV 
source which occurs during the thermahzation stage of 
the reaction. At mid-rapidity there is Httle contribution 
from the PLF source remaining. There is, however some 
observed contribution from the TLF source at low veloc- 
ities in the IV source frame. In the following analysis 
yields assigned to the TLF source have been subtracted 
from the experimental yields. 

As in reference > the AMD- V model calculations 0, 
0| indicate a significant slowing in the rate of kinetic en- 
ergy change near a velocity of 3.5 cm/ns. This signals 
the end of the IV (or pre-equilibrium) emission stages 
and entry into the region of slower nuclear de-excitation 
modes, i.e., evaporation, fission and/or fragmentation. 
At that point the sensitivity of the emission energy to 
time is significantly reduced. Therefore we take the tem- 
perature at the time corresponding to the velocity of 3.5 
cm/ns to be that of the hot nucleus at the beginning of 
the final statistical emission stage. 

We present, in Figure 1, the double isotope ratio tem- 
peratures as a function of emission time. As time in- 
creases, each of the temperature evolution curves rises 
to a maximum and then decreases. Maximum temper- 
atures of 10-15 MeV are observed at times in the range 
of 95 to 110 fm/c. After that the temperatures drop 
monotonically. Since the transport model calculations 
indicate that contributions from the late stage evapora- 
tion become dominant below Ysurf = 3.5 cm/ns (see 
also reference TtJj) we terminate the curves at the times 
corresponding to that velocity. We take these times to 
represent the end of the IV (or pre-equilibrium) emission 
stage. It is worth noting that the corresponding temper- 
atures, presented in column 4 of Table 1, are near 6 MeV 
and thus very similar to the limiting temperatures pre- 
viously derived from a systematic investigation of caloric 
curve measurements [22(, in this mass region. 

The trends in Fi gure 1 are very similar to those re- 
ported in reference \MMj for the temperature evolution 
in the reactions of 26-47 MeV/u ^^Zn projectiles with 
Ni, Mo and Au targets. The possible interpretation of 
these trends is extensively discussed in that reference. It 
is concluded that the data are consistent with chemical 
and thermal equilibration of the ensemble of events sam- 
pled, at least on a local basis. 

If thermal and chemical equilibrium are achieved by 
the time of entry into the final fragmentation or evapora- 
tion phase of the reaction and the density at that point is 
not too high [2^ an agreement between the thermal tem- 
perature and the double isotope ratio temperature Thhc 
at that time can be expected. The TLF source fit slope 
parameters, Tsiope for alpha particles are seen in column 
5 of Table 1 to be slightly lower than the latest time Thhc 
temperatures. This is not surprising since these Thhc 
temperatures represent the temperature at the start of 
the de-excitation cascade while the TLF source fits to the 
alpha spectra can be expected to return only an apparent 
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FIG. 1: Time evolution of the double isotope ratio tempera- 
ture, Thhe, for the four different reactions. 35 MeV/u ^*Zn-|- 
^^Mo- red diamonds, 40 MeV/u ''°Ar -I- "^Sn- orange trian- 
gles, 55 MeV/u ^'^Al + ^^^Sn- yellow circles and 47 MeV/u 
^^Zn + ''^Mo - blue squares. 



temperature reflecting the entire evaporation and coo ling 
history of the TLF source. In previous work [2^ |2y, [23 
we have found that the slope parameters for the alpha 
particle emission from the TLF source most closely ap- 
proximates the intial thermal temperature of this source, 
reflecting the higher fraction of the alpha emission in the 
earlier part of the de-excitation cascade. The agreement 
between these thermal fit apparent temperatures and the 
late time chemical temperatures obtained from the iso- 
tope ratio temperatures, as presented in column 4 is quite 
reasonable given the associated uncertainties. 



DE-EXCITATION CASCADES 

If, in fact, the hot nuclei have achieved equilibration 
at very similar excitation energies and temperatures, this 
should be reflected in their de-excitation patterns. Figure 
2 presents the hot nucleus masses, excitation energies and 
Thhc temperatures for the four different systems (in row 
1 of the Figure) and a comparison of average slope tem- 
peratures, multiplicities and kinetic energies for n, p, d, t, 
■^He and *He emission during the TLF de-excitation. For 
each of the quantities presented in the figure, a dashed 
line representing the average values obtained for the 35 
MeV/u 64Zn-f ^^Mo, 40 MeV/u ^^Ar + "^gn and 47 
MeV/u ^''Zn + ^^Mo reactions ( which lead to nuclei 
with similar N/Z ratios) is also shown. For each individ- 
ual system, deviations from those averages are seen to be 
relatively small except for the enhanced neutron multi- 
plicity for the 55 MeV/u ^^Al + ^^^Sn reaction, which 
is to be expected given the larger total excitation energy 
and significantly higher N/Z of the hot nucleus produced 
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Reaction 


Acomp 


E* 
MeV/u 


MeV 


slope 

MeV 




E7AA 
MeV/u 


35 MeV/u ^^Zn+^^Mo 
35 MeV/u ""Ar+^^Sn 
47 MeV/u ^''Zn+s^jyto 

55 MeV/u 2''Al+^2'*Sn 


96.4 ±8.7 
107±9.6 

85.4 ±7.8 
114 ± 10 


4.88 ± 0.49 
4.98 ±0.50 
5.14 ±0.51 
5.31 ±0.53 


5.93 ±0.59 
5.98 ±0.60 
5.60 ±0.56 
6.00 ±0.60 


5.10 ±0.50 
5.80 ±0.30 
5.50 ±0.03 
5.60 ±0.10 


38.7 ±3.5 
41.0 ±3.7 
35.5 ±3.3 
44.5 ±4.1 


12.1 ± 1.7 
13.0 ±1.8 
12.4 ± 1.8 
13.6 ±1.9 



TABLE I; Properties of the Hot Nuclei Produced in mid-central collisions. Errors are estimated absolute uncertainties. 



with this neutron rich target [l^ . Since the total excita- 
tion energies do vary, we present, in column 6 of Table 1, 
the total mass removed by TLF ejectiles and in column 
7 , the ratio of total excitation energy to total mass re- 
moved by the TLF ejectiles. This latter number is seen 
to be quite stable within experimental uncertainties. The 
average value, 12.8 ± 1.8 MeV, is close to values obtained 
in an early study of the energetics and de-excitation in 
reactions of 27 MeV/u ^°Ar with Ag and Ho leading to 
somewhat heavier hot nuclei with similar excitation en- 
ergies In that study, which detected heavy residues, 
equilibration and isotropy of emission was inferred from 
the residue velocity and angular distributions. In the 
present work, isotropy of emission from the TLF source is 
indicated by the success of that assumption in the source 
fitting process. Except for the expected additional favor- 
ing of neutron emission over charged particle emission 
for the reaction with the ^^^Sn target, the very similar 
de-excitation patterns and energetics for the four hot nu- 
clei which were produced in somewhat different dynam- 
ical conditions, is most easily understood as a reflection 
of statistical emission from very similar equilibrated hot 
nuclei. 

It is interesting therefore to ask how the values in Fig- 
ure 2 compare with those calculated using current statis- 
tical model codes. For this purpose we have carried out 
calculations using both the GEMINI code |2^ and the 
Statistical Multi-fragmentation Model code (SMM) 
Except for an exploration of the effect of varying the level 
density parameter, a, these codes were used with normal 
default parameters and no attempt has been made to 
tailor them to the current data. For these highly excited 
systems the change in the multiplicities and average en- 
ergies which results from varying the GEMINI code level 
density parameter was found to be rather small. The 
GEMINI results in Figure 2 were obtained with a level 
density parameter, a = A/9. Yoshida et al. found this to 
be a reasonable average value in their study of neutron 
emission in the reactions of ^°Ar ± ^■^Mo (l2l | . 

As is seen in that figure, the calculated relative val- 
ues of the multiplicities and average energies for each 
system show trends quite similar to those of the exper- 
iments, suggesting that the statistical picture has some 
validity. However, the statistical model calculations sys- 
tematically produce more Z=l ejectiles and fewer neu- 
trons than have been observed in the experiments. A 



similar under-prediction of neutron emission is seen in 
reference In addition, except for the protons, the 

calculated average kinetic energies of the charged ejec- 
tiles are somewhat higher than the experimental values. 
Given that the dynamic calculations indicate that the 
initial hot de-exciting nuclei produced in these collisions 
have more complicated shapes and density distributions 
than those assumed in the calculations it is perhaps not 
so surprising that some differences between the experi- 
mental results and the calculated results are observed. 



CONCLUSION 

Establishing that equilibration has occurred in Fermi- 
energy collisions is difficult, but remains an important 
pre-requisite in determining the applicability of many of 
the techniques which have been utilized or proposed as 
tools to derive information on the nuclear equation of 
state at non normal densities 0, 0, 0- ^^^^ 
per we have presented results of a Ghoshal-like exper- 
imental test of equilibration of medium mass hot com- 
posite nuclei produced in central collisions of 35 to 55 
MeV/u projectiles with medium mass targets. Although, 
in contrast to the low energy experiments of Ghoshal, the 
present entrance channel mass asymmetries and energies 
differ significantly for the four different cases studied and 
the masses of the final residues are less than half of the 
entrance channel masses, comparisons of the early tem- 
perature evolution and of the subsequent de-excitation 
patterns for the four systems indicate that very similar 
hot nuclei are indeed produced and that they decay in a 
very similar manner. Multiplicities and energies derived 
from statistical model calculations of the de-excitation 
of these hot nuclei show system to system trends simi- 
lar to those of the experimental results but some quan- 
titative differences in the competition between neutron 
and light charged particle emission. Overall, the results 
are most easily understood as a reflection of statistical 
emission from very similar equilibrated hot nuclei. That 
such systems achieve equilibration is in agreement with 
the conclusions of a recent study in which evidence for 
equilibration was derived from observations of the tem- 
perature evolution in similar collisions |lOl |. 
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FIG. 2; Comparison of initial properties of the hot nuclei and de-excitation modes for the four different systems, (a) Masses of 
the de-exciting hot nuclei, (b) Excitation Energies of the hot nuclei, (c) Double Isotope Temperatures, Thhb, (d) - (f ) Slope 
temperatures for n,p and , (g) - (1) multiplicities of n, p, d, t, ^He and , (m) - (r) average kinetic energies of n, p, d, t, ^He 
and in the hot nucleus frame. Solid points represent the data for the four different systems. Open squares show the results of 
Gemini calculations (22|. Open triangles show the results of SMM calculations [s^. 



ACKNOWLEDGEMENTS 



This work was supported by the United States Depart- 
ment of Energy under Grant # DE-FG03- 93ER40773 
and by The Robert A. Welch Foundation under Grant # 
A0330. 



* E-mail at wangQcomp.tamu.edu 

t Now at FNRS and IPN, Universite Louvain-Neuve, Bel- 
gium Catholique de Louvain, B-1348 

* on leave from Shanghai Institute of Nuclear Research, 
Chinese Academy of Sciences, Shanghai 201800, China 

^ Now at University of Kansas, Lawrence, Kansas 66045- 
7582 

^ Now at LCP Caen, ISMRA, IN2P3-CNRS, F-14050 



7 



Caen, France 
[1] S. N. Ghoshal, Phys. Rev. 80, 939 (1950). 
[2] N. Bohr, Nature 137, 344 (1936). 

[3] J. Blatt and V. F. Weisskopf, Theoretical Nuclear Physics 

(John Wiley & Sons, Inc., New York, 1952). 
[4] E. Suraud, C. Gregoire and B. Tamain, Prog. Part, and 

Nucl. Phys. 23, 357 (1989). 
[5] B. Tamain and D. Durand University of Caen Report 

LPCC 96-16 (1996) and references therein. 
[6] P. Chomaz, To appear in Proceedings of the INPC 2001 

Conference, Berkeley, Ca., July 2001 
[7] J.B. Elliott et al, Phys. Rev. C67 024609 (2003). 
[8] A. Ono, Phys. Rev. C 59, 853 (1999) 
[9] R. Wada, T. Keutgen, K. Hagel, Y. G. Ma, J. Wang, 

M. Murray, L. Qin, P. Smith, J. B. Natowitz, R. Al- 

farro, J. Cibor, M. Cinausero, Y. El Masri, D. Fabris, E. 

Fioretto, A. Keksis, M. Lunardon, A. Makeev, N. Marie, 

E. Martin, A. Martinez-Davalos, A. Menchaca-Rocha, G. 
Nebbia, G. Prete, V. Rizzi, A. Ruangma, D. V. Shetty, 
G. Souliotis, P. Staszel, M. Veselsky, G. Viesti, E. M. 
Winchester, S. J. Yennello, Z. Majka and A. Ono, Phys. 
Rev. C 69 044610 (2004). 

[10] J. Wang, R. Wada, T. Keutgen, K. Hagel, Y. G. Ma, 
M. Murray, L. Qin, A. Botvina, S. Kowalski, T. Ma- 
terna, J. B. Natowitz, R. Alfarro, J. Cibor, M. Cinau- 
sero, Y. El Masri, D. Fabris, E. Fioretto, A. Keksis, M. 
Lunardon, A. Makeev, N. Marie, E. Martin, Z. Majka, 
A. Martinez-Davalos, A. Menchaca-Rocha, G. Nebbia, 
G. Prete, V. Rizzi, A. Ruangma, D. V. Shetty, G. Souli- 
otis, P. Staszel, M. Veselsky, G. Viesti, E. M. Winch- 
ester, S. J. Yennello, W. Zipper, A. Ono, ArXiV Preprint 
nucl-ex/0408002 (2004) Submitted to Physical Review C. 

[11] L Tilquin, Y. El Masri, M. Parlog, P. CoUon, M. Hadri, 
T. Keutgen, J. Lehmann, P. Leleux, P. Lipnik, A. Ninane, 

F. Hanappe, G. Bizard, D. Durand, P. Mosrin, J. Peter, 
R. Rgimbart and B. Tamain, Nucl. Instr. Meth. A 365 
446 (1995). 

[12] K. Yoshida, J. Kasagi, H. Hama, M. Sakurai, M. Ko- 
dama, K. Furutaka, K. leki, W. Galster, T. Kubo, M. 
Ishihara, and A. Galonsky, Phys. Rev. C 46, 961-971 
(1992) 

[13] J. Cibor, R. Wada, K. Hagel, M. Lunardon, N. Marie, R. 

Alfaro, W. Q. Shen, B. Xiao, Y. Zhao, J. Li, B. A. Li, M. 

Murray, J. B. Natowitz, Z. Majka and P. Staszel, Phys. 

Lett. B 473, 29 (2000). 
[14] K. Hagel, R. Wada, J. Cibor, M. Lunardon, N. Marie, R. 

Alfaro, W. Shen, B. Xiao, Y. Zhao, Z. Majka, J. Li, P. 

Staszel, B.-A. Li, M. Murray, T. Keutgen, A. Bonasera, 

and J. B. Natowitz, Phys. Rev. C 62, 034607 (2000). 
[15] J. Cibor et al., in "Isospin Physics in Heavy-Ion Collisions 

at Intermediate Energies", Eds.Bao-An Li and W. Udo 

Schroeder, NOVA Science Publishers, Inc. (New York), 

(2001) 

[16] Y. G. Ma, J. B. Natowitz, R. Wada, K. Hagel, J. Wang, 



T. Keutgen, Z. Majka, M. Murray, L. Qin, P. Smith, R. 
Alfaro, J. Cibor, M. Cinausero, Y. El Masri, D. Fabris, E. 
Fioretto, A. Keksis, M. Lunardon, A. Makeev, N. Marie, 
E. Martin, A. Martinez-Davalos, A. Menchaca-Rocha, G. 
Nebbia, G. Prete, V. Rizzi, A. Ruangma, D. V. Shetty, 
G. Souliotis, P. Staszel, M. Veselsky, G. Viesti, E. M. 
Winchester, S. J. Yennello, nucl-ex/0410018 Submitted 
to Physical Review C. 
[17] T.C. Awes, G. Poggi, C.K. Gelbke, B.B. Back, B.C. 
Glagola, H. Breuer and V.E. Viola, Jr., Phys. Rev. C 
24, 89 (1981). 

[18] T. C. Awes, G. Poggi, S. Saini, C. K. Gelbke, R. Legrain 
and G. D. WestfaU, Phys. Lett. 103B, 417 (1981) 

[19] D. Prindle, A. Elmaani, C. Hyde- Wright, W. Jiang, A. 
A. Sonzogni, R. Vandenbosch, D. Bowman, G. Cron, P. 
Danielewicz, J. Dinius, W. Hsi, W.G. Lynch, C. Montoya, 
G. Peaslee, C. Schwarz, M. B. Tsang, C. Wilhams, R.T. 
de Souza, D. Fox, and T. Moore, Phys. Rev. C 57, 1305 
(1998). 

[20] R. Wada, D. Fabris, K. Hagel, G. Nebbia, Y. Lou, M. Go- 
nin, J. B. Natowitz, R. Billerey, B. Cheynis, A. Demeyer, 
D. Drain, D. Guinet, C. Pastor, L. Vagneron,K. Zaid, J. 
Alarja, A. Giorni, D. Heuer, C. Morand, B. Viano, C. 
Mazur, C. Ngo, S. Leray, R. Lucas, M. Ribrag and E. 
Tomasi, Phys. Rev. C 39, 497 (1989) 

[21] K. Siimmerer and B. Blank, Phys. Rev. C 61, 034607 
(2000) 

[22] S. Kowalski et al., work in progress. 

[23] J.B. Natowitz, R. Wada, K. Hagel, T Keutgen, M. Mur- 
ray, Y. G. Ma , A. Makeev, L. Qin , P. Smith and 
C.Hamihon, Phys. Rev. C 65, 034618 (2002). 

[24] L. G. Sobotka, R. J. Charity, J. Toke, and W. U. Schroder 
Phys. Rev. Lett. 93, 132702 (2004). 

[25] Z. Majka, P. Staszel, J. Cibor, J. Natowitz, K. Hagel, J. 
Li, N. Mdeiwayeh, R. Wada and Y. Zhou, Phys. Rev. C 
55, 2291 (1997). 

[26] K. Hagel, D. Fabris, P. Gonthier, H. Ho, Y. Lou, Z. 
Majka, G. Mouchaty, M.N. Namboodiri, J.B. Natowitz, 
G. Nebbia, R.P. Schmitt, G. Viesti, R. Wada and B. 
Wilkins, Nucl. Phys. A486, 429 (1988). 

[27] M. Gonin, L. Cooke, K. Hagel, Y. Lou, J. B. Natowitz, R. 
P. Schmitt, S. Shlomo, B. Srivastava, W. Turmel, H. Ut- 
sunomiya, R. Wada, G. Nardelli, G. Nebbia, G. Viesti, R. 
Zanon, B. Fornal, G. Prete, K. Niita, S. Hannuschke, P. 
Gonthier, and B. Wilkins, Phys. Rev. C 42, 2125 (1990). 

[28] M.F. Rivet, B. Borderie, H. Gauvin, D. Gardes, C. 
Cabot, F. Hannape and J. Peter, Phys. Rev. C 34, 1282 
(1986). 

[29] R. J. Charity, et al. Nucl. Phys. A483, 371 (1988) 
http://www.chemistry.wustl.edu/faculty/charity.html 

[30] J.P. Bondorf, A.S. Botvina, A.S. Iljinov, I.N. Mishustin, 
K. Sneppen. Phys. Rept. 257 133 (1995). A. Botvina, 
private communication. 



